Pump process of the rotatory molecular motor and its energy efficiency 
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The pump process of the ratchet model inspired by the Fo rotatory motor of ATP synthase 
is investigated. In this model there are two kinds of characteristic time. One is dynamical, the 
relaxation time of the system. Others are chemical, the chemical reaction rates at which a proton 
binds to or dissociates from the motor protein. The inequalities between them affect the behavior 
of the physical quantities, such as the rotation velocity and the proton pumping rates across the 
membrane. The energy transduction efficiency is calculated and the condition under which the 
efflciency can become higher is discussed. The proton pumping rate and the efficiency have a peak 
where a certain set of inequalities between the chemical reaction rates and the reciprocal of the 
relaxation time holds. The efficiency also has a peak for a certain value of the load. The best 
efficiency condition for the pump process is consistent with that for the motor process. 

PACS numbers: 05.40.-a, 05.60.-k, 87.15.-v 



I. INTRODUCTION 

Biomolecular motors are kinds of protein molecules 
which generate a unidirectional motion in the situation 
that many surrounding molecules, e.g. water molecules, 
collide with the protein molecule many times due to ther- 
mal fluctuation They play many essential roles in 
life; muscle contraction, transport in cells, etc. It is inter- 
esting and important question how motor proteins gen- 
erate a unidirectional motion. It is natural to think that 
they utilize this random fluctuation by the collision, al- 
though there has been no unambiguous evidence reported 
so far that this is essential. 

Ratchet systems, which take out a net finite current 
from random and undirected noise, have been intensively 
investigated for modeling molecular motors As known 
by the famous Feynman's lecture^ and many investiga- 
tions thereafter!^ 0,0, Hi, a net finite directed motion is 
the consequence of the violation of the symmetry of the 
system, e.g. the spatial form of the potential. In mod- 
eling molecular motors, the system also needs to violate 
the detailed balance between several different states, i.e. 
conformations. The detailed balance is violated by the 
undirected chemical energy supply, e.g. the ATP hy- 
drolysis. Ratchet systems are applied to various physical 
systems other than modeling molecular motors ja 0, 

The Fo part of the ATP synthase is known as a rota- 
tory molecular energy transducer Hol ITll IT^ . It is 
embedded in the membrane, and transduces the energy of 
the transmembrane proton concentration gradient, when 
a proton passes through the membrane via the Fq part, 
into that of the rotation of its internal subunits in the 
motor process. The ratchet model of it was presented 
qualitatively by Junge, Lill and Engelbrecht[I3- Elston, 
Wang, and Oster improved this model by taking account 
of the electrostatic interaction between the residues Ts^ . 
This ratchet model has an advantage that the coupling 
between the rotation and the proton flow is easy to see, 
therefore the energy transduction efficiency is easily esti- 



mated. In this model, there are two important kinds of 
characteristic time, dynamical and chemical. The former 
is the relaxation time of the convection-diffusion system. 
The latter are the chemical reaction rates at which a pro- 
ton binds to or dissociate from the motor protein when 
it passes through the membrane. We investigated the ef- 
ficiency in detail and pointed out that in the motor pro- 
cess, the inequalities between them affect the behavior of 
it [l^- For the motor to work most efficiently, certain 
inequalities between them need to hold. 

The Fo part also works as a pump^3|; it gains en- 
ergy from the (inverse) rotation and pumps out protons 
against the transmembrane proton gradient. Which pro- 
cess this protein adopts, motor or pump, under a given 
condition and what controls this determination are im- 
portant questions, but satisfying answers to these have 
not been obtained so far. 

In this paper we investigate the physical quantities in 
the pump process, such as the rotation velocity, the pro- 
ton pumping rate, and the energy transduction efficiency 
The condition under which the protein works as a pump 
with higher efficiency is discussed. We will compare the 
most efficient condition for the pump process with that 
for the motor process. 



II. MODEL 

The ATP synthase is schematically shown in FIG^a). 
The upper side of the figure is the basic (inner) side of 
the cell and the lower side is the acidic (outer) side. The 
proton concentration of the acidic side is kept higher than 
that of the basic side by respiratory chains. The Fo part 
is composed of the subunits denoted by the Roman let- 
ters, a, b, and c, and the Fi part by the Greek letters, 
a, (3, 7, and 5. The rotating ring (the c-ring)is com- 
posed of the assembly of the c-subunits, each of which 
has a proton binding site, carboxyrate Asp61, almost in 
the middle of itself. The a-subunit is fixed in the mem- 
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brane. The proton channel is the interface between the 
a-subunit and the c-ring. There are two c-subunits and 
two proton paths in the channel. One of the paths con- 
nects the left binding site of the c-subunits to the ba- 
sic side and the other connects the right binding site to 
the acidic side. The other c-subunits are in the mem- 
brane. The proton binding site can be both protonated 
and unprotonated only if it is in the channel. Otherwise it 
can be only protonated due to the hydrophobicity of the 
membrane. The conformational change of the c-subunits 
between the two states, protonated and unprotonated, is 
small ^-iid neglect the effect on the place of the 
binding sites. The states of the c-ring are determined 
by protonations of the two proton binding sites in the 
channel; the empty(E) state with the both binding sites 
unprotonated, the right(R) state with the right binding 
site protonated and the left binding site unprotonated, 
the left(L) state with the left binding site protonated 
and the right binding site unprotonated, and the full(F) 
state with the both binding sites protonated. As a motor 
the c-ring rotates rightward and as a pump leftward in 
FIGH The 7- subunit connects the Fo part to the Fi 
part. 

We investigate this system by using the 'simply bi- 
ased diffusion model' ITEL It works as follows: 
The c-ring moves randomly by collision with surround- 
ing molecules and is forced to rotate by the external load 
from the Fi part. But the hydrophobic interaction pre- 
vents the unprotonated binding site from moving into the 
membrane, i.e. the E and R states cannot move leftward 
out of the channel and the E and L states not rightward. 
If the chemical reaction rates, at which a proton binds 
to or dissociates from the binding sites in the channel, 
are tuned appropriately, it becomes possible to take out 
a net motion against the load or to pump protons against 
the concentration gradient. 

We consider the pump process of this system. A proton 
passes through the channel from the basic side to the 
acidic side as follows: First a proton in the basic side 
flows into the channel and binds to the left binding site. 
Next the proton goes through the membrane with the 
leftward rotation of the c-ring. Finally the proton comes 
back to the channel, dissociates from the right binding 
site, and flows out to the acidic side. In this way, a proton 
is pumped out by the rotation of the c-ring driven by 
the external load, which is gained by the energy of the 
ATP hydrolysis in the external part. We assume that 
in the pump process, the external load exerted is larger 
than that in the motor process, so that the c-ring can 
rotate leftward. In this sense, the external 'load' should 
be exactly called 'driving force' of the rotation in the 
pump process. However, we will use the term 'load' in 
the pump process since we deal with the same model as 
in the motor process. 

First of all, we give a qualitative consideration of the 
energy transduction, the coupling between the rotation 
and the proton flow. The good process in the channel 
is shown in FIGI21 A proton binds to the left binding 
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FIG. 1: The upper figure (a) is the ATP synthase. When it 
operates, the shaded parts (7, e, and the c-ring) rotate. In 
the pump process, the external load force from Fi rotates the 
c-ring in the direction of the arrow and pumps out a proton 
from the upper to the lower side. This direction of rotation 
is opposite to that in the motor process. In the lower figures, 
(b), the states of the c-ring are shown. The rectangle frame is 
the proton channel viewed from the a-subunit and the shaded 
region is the c-ring. The states depend on the protonations 
of the two proton binding sites in the channel. The figure (a) 
corresponds to the R-state. 



site in the left figure. Then the c-ring moves one step 
leftward in the center figure. Finally the proton bound 
to the right binding site dissociates in the right figure 
and the process returns back to the start. In this process 
one step of the c-ring rotation corresponds to one proton 
pumping. This is the efficient energy transduction. A 
bad process is shown in FIG|21 In the upper-right fig- 
ure, a proton binds to the left binding site before the 
proton bound to the right binding site dissociates. Then 
in the lower-right figure, the c-ring moves leftward with- 
out proton dissociation. Therefore more than one step of 
the rotation corresponds to one proton pumping. In this 
process the energy is transduced inefficiently. Another 
bad process is shown in FIG0] In the upper-middle fig- 
ure, a proton binds to the right binding site from the 
acidic side before the rotation of the c-ring. Next the 
c-ring rotates leftward in the upper-right figure. Finally 
the proton bound to the left binding site dissociates then 
the c-ring rotates rightward due to diffusion. In this pro- 
cess, a proton passes through the membrane along the 
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:rotation by load 



rotation by load 



FIG. 2: The good process of the coupling between the proton 
pumping and the motion of the c-ring. In this process, one 
step of the rotation tightly couples to one proton pumping. 
The energy is transduced efficiently. 




•< — *• :rotation by load or diffusion 

FIG. 3: A bad process. One proton pumping corresponds to 
more than one step of the rotation. The energy of the motion 
is wasted in this process. 



FIG. 5: The Fo part viewed from the basic side. The position 
variable x is defined as a rotation angle of the c-ring. 



PE{x,t) 
PF{x,t) 



(3) 



FIG. 4: Another bad process. A proton flows in the undesir- 
able direction without rotation. This process is possible since 
the proton concentration of the outer side is kept higher than 
that of the inner side. This process suppresses the efficient 
energy transduction. 



gradient, in the undesirable direction, without net step 
of the c-ring. These bad processes generate the loss of 
the energy transduction. For the higher efficiency, the 
good process should become more and the bad processes 
should become less. 

The system is described by the Fokker-Planck equa- 
tion. 



|p(x,i) = -^n(^,i) 



nix,t) ^ -f{{x)p{x,t) 



-K{x) ■pix,t), 



i^^p(x,t). 



(1) 



(2) 
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nF{x,t) 



(4) 



Here p{x, t) is 4-component probability where Pi{x, t) de- 
scribes the probability that the c-ring in state i is at posi- 
tion x at time t. The state index i refers to the state men- 
tioned before; E, R, L, and F. Similarly, 11 describes the 
flow of the probability. K is the transition matrix which 
describes changes between the states and will be given 
later. D and 7 are the diffusion constant and the friction 
constant, respectively. They satisfy the Einstein's rela- 
tion, D = jk^T where T is the temperature and is 
the Boltzmann constant. The matrix f = diag[r, r, r, r] 
is the external load. The position variable x is defined 
as the rotation angle with the origin placed at the cen- 
ter of the c-ring and takes its value from to S = 2tt/N 
where N denotes the number of c-subunit consisting in 
the c-ring. (See FigEl) Hereafter, we fix N =12, which 
is the case of E. Coli. The position a; = is defined as the 
situation where the left binding site is at the left bound- 
ary of the channel and a; = (5 is as that where the right 
binding site is at the right boundary of the channel. 

As mentioned earlier, the conformational change of the 
c-subunit due to the protonation of the binding site is 
small. Therefore we have assumed that all the state take 
the common values of the diffusion constant D and the 



4 



friction constant 7. As an effect of the conformational 
change, we wih only take account of the hydrophobic- 
ity, which determines whether or not a c-subunit in the 
channel can move into the membrane. 

The system has a periodicity in the sense that when 
one binding site in the channel goes out and moves into 
the membrane, another binding site comes into the chan- 
nel. This periodicity and the hydrophobicity impose the 
boundary conditions on the flows. 



PFiS)=PFiO). 



The transition matrix K is given as 



K = 



(10) 



and, 



nE(o,i) = nE(<5,t) = 0, 



nL(o,t) = nR,(5,i), 



(5) 



(6) 



(7) 



(8) 



Similarly, the following boundary conditions of the prob- 
abilities are imposed. 
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(11) 

where fc*(i=R,L,j=in,out) denotes the rate at which a 
proton binds to (in) or dissociate from (out) the right (R) 
or left(L) binding site. Here we assumed that (1) the 
chemical reaction is sufficiently faster than the motion of 
the c-ring; (2) there is no proton hopping between the 
binding sites and no correlation between the reactions; 
(3) the reaction rates are independent of the position, x. 

The steady state solution is of interest therefore the 
time variable is omitted hereafter. Under these assump- 
tions, this system can be solved analytically. The matrix 
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diagonalizes the transition matrix K 



Q ^KQ = 



(^in + ^out) 



Then Eq.Q in the steady state is reduced to 
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(Q-ip(x)). (14) 



The solution is 



(Q-ip(a:)) = 



where (i = 1, 2, • • • , 8) are integral constants and 



(15) 
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TABLE I: Parameters used for calculation. 



^_ 7r±V(7r)^+4ZJ(fc^^ + Cty 



'yLR - 2Z? 

(19) 

The integral constants C^'s are determined by the bound- 
ary conditions of the probabilities and their flows, 
Eas.lf5|)- (|10|l . To close the algebraic equations for C^'s, 
the normalization condition of the probability 



(20) 



is needed. 

Physical quantities which we will investigate are calcu- 
lated as follows: First, the average rotation velocity {v) 
is given as, 

(«) = En^' (21) 

i 

dxIVi{x) 

^Tp, + D[pm~P^m, (22) 
dxpi(x). (23) 



where 



and. 



Next, the proton pumping rate A^p(H+) is written as, 



^p(H+) = -(Jr + Jl) 



(24) 



where, 



= -koutPB. + k^,pE - k^^^pF + k^pL, (25) 



(26) 



Here Jr and Jl describe the proton flow of the right and 
left binding site, respectively. Finalh^ let us define the 
efficiency of the energy transduction 19]. The energy in- 
put per unit time is defined as the product of the external 
load torque r and the rotation velocity, (v). The energy 
output per unit time is the product of the free energy 
due to the transmembrane proton gradient, AG, and the 
proton pumping rate, A^p(H+). Therefore the efficiency 
e is defined as, 



_ AG ■ jVp(H+) 

T-{v) 



(27) 
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III. RESULTS 

We investigate the dependence of the quantities, (w), 
iVp(H+), and e on the transition rates, K. The transition 
rates are given by 



'^in 
"■in 
"-out 
'^out 



= 10 



K 



IQ-pHBg't'/kBT 
lQ-pK^^-V/2kBT 
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(28) 



The parameters used for the calculation are, according to 
Ref.na, given in TABLE^ They areplausible values in 
living things. In our previous work |lfil |. we pointed out 
that the inequalities between the chemical reaction rates, 
k'j 's, each other and the reciprocal of the relaxation time 
of the system, T^^^^ = {'^tY /D, affect the behavior of 
the physical quantities in the motor process. It is natural 
to expect that these inequalities also play important roles 
in the pump process. The chemical reaction rates depend 
on many kinds of factor, e.g. the diffusion constant of 
the proton, the acidity of the proton binding site, the 
concentrations of proton of both sides, etc. The overall 
factor IQ-'^ means the effective proton absorption rate 
of the path, which is regarded as common for all the 
transition rates. 

The free energy accompanying the proton transloca- 
tion due to the transmembrane proton gradient is given 
as 

AG = V + k^T\n[lO^P"] 



31.6[pN • nm] 



(29) 



where ApiJ = pHs — pHa- The detailed balance is vio- 
lated by this free energy and the pumping is against this 
energy. 

We investigate in detail two sets of inequalities between 
the transition rates; case 1) /c^^j > fc^^^ > kj^ > fc?^ 
and case 2) k^^^^ > kj^ > fc^^ > kj"^. They are derived 
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FIG. 6: The rotation velocity {v) in the case of the set of 
inequahties fc^ut > k^ut > > ^li {pKa—4.5, case 1)). 
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FIG. 7: The rotation velocity (u) in the case of the set of 
inequalities k^ut > kj^ > k^^^t > fcii (^^^0=5. 5, case 2)). 



from the parameters in TABLE U and the condition un- 
der which our model works well as a motor. We take the 
value of the acidity of proton binding site, pKa=4:.5 and 
5.5 as typical value which gives the set of inequalities 1) 
and 2), respectively. 

The external load t is also regarded as a controllable 
parameter, since the load has not been so far measured 
experimentally in the pump process we discuss here. 

The rotation velocity (v) is plotted in FIGs|Hland[71in 
the case 1) and 2), respectively. There is no qualitative 
difference between two cases. There is almost no rota- 
tion when K is small, i. e. the chemical reaction rates are 
low. Then it monotonically increases and finally satu- 
rates at K >15, which is out of the figures. As expected, 
it becomes faster for the larger load. 

The proton pumping rate iVp(H+) is plotted in FIGsEI 
and 1^1 Almost no pumping is observed for small K in 
both cases. Then it increases as K increases. It reaches 
a peak at a certain finite value of K, such that fc^^^ ~ 
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FIG. 8: The proton pumping rate A^p(H"'") in the case 1). For 
fixed load r, this quantity has a peak at about k^^^ ~ = 
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FIG. 9: Afp(H+) in the case 2). 



^reiax' ^^"^ then falls down abruptly. The larger the load, 
the more the peak value increases. 

The energy transduction efficiency e is plotted in 
FIGs^l and ^2 For respective fixed loads r, it has 
a peak for a certain finite value of K. This fact re- 
flects the above behaviors of the rotation velocity and 
the proton pumping rate. The peak is placed on K such 



that T.-;, > > Tj^^, 



^ holds. Here Tioad = 
is the characteristic time that it takes for the c-ring to 
rotate one step leftward by the load. The variation of 
the peak value with respect to the load r is shown in 
FIGs^l and ^1 For the variation of the load r, there 
exists an optimal value of r such that the peak value 
of the efficiency becomes maximum. In the case 1), 
T ~75pN-nm and r ^90pN-nm in the case 2). Note that 
the efficiency in the case 1) is about one order larger 
than that in the case 2). Therefore, the pump mecha- 
nism works most efficiently when the set of inequalities 

^ fcout > T-L > ^out > k^n > holds and the 
optimal value of t ~75pN-nm is fulfilled. 

We can understand qualitatively the reason why the 
above inequalities are required for the efhcient pump pro- 
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FIG. 10: The energy transduction efficiency e in tlie case 1). 
For a fixed load r, there is a pealt at about T~J^^^ > k^^t > 

^ load • 
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FIG. 11: The energy transduction efficiency e in the case 2). 
For a fixed load r, there is a pealt at about T~J^^^ > k]^^^ > 
^load- this case the value is about one order smaller than 
the former case. 



cess. Firstly, the inequality T~^^^ > k^^^ is required to 
suppress the bad process shown in FIG01 If the inequal- 
ity is not satisfied, the left binding site in the channel 
becomes unprotonated before the relaxation by the fluc- 
tuation hence it cannot move into the membrane. As a 
result the bad process in FIG^ can occur easily. Sec- 
ondly, the inequality /c^^^ > T^^^^ indicates that the bad 
process shown in FIGOl rotation without proton pump- 
ing, rarely occurs. Thirdly, the inequality fc^^^ > k^^, 
which distinguishes the case 1) from 2), is needed to pre- 
vent the c-ring from moving rightward. 

For Topt, we note that it is determined in highly non- 
trivial manners. To see this, let us go back to Ea. l27|) . As 
we see in FIGs l6l9l both the numerator of Ea. H27|) . which 
is proportional to the proton pumping rate iVp(H+), and 
the denominator, which is proportional to the rotation 
velocity (v) and the load r, increase monotonically as the 
load increases. Therefore, the subtle difference between 
the increasing rates of them determines Topt- 
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FIG. 12: The variation of the peak value of the energy trans- 
duction efficiency e in the case 1) with respect to the load 
r. Below the critical load Tc ~60(pN-nm), the system cannot 
work as a pump. It has a peak for a certain finite value of r, 
Topt ~75(pN-nm). 
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FIG. 13: The variation of the peak value of the energy trans- 
duction efficiency e in the case 2) with respect to the load 
r. For r < Tc ~70(pN-nm), the system cannot work as 
a pump. There is a peak for a certain finite value of r, 
Topt ~90(pN-nm). The value is about one order smaller than 
the case 1). 



IV. SUMMARY AND DISCUSSION 

We have investigated the pump process of the ratchet 
model of the rotatory molecular motor. The model can 
be solved analytically and the physical quantities such as 
the rotation velocity and the proton pumping rate are cal- 
culated. The energy transduction efficiency is obtained 
from these physical quantities. The best efficiency condi- 
tion is that the set of inequalities > k^^^ > T{^1^ > 
^out > ^hi > ^in holds undcr an optimal value of the load 

''opt- 

We compare the result obtained here with that for the 
motor process. In Ref. 16] we find that the most effi- 
cient condition for the motor process is that the set of 
inequalities k^^, > T~l^ > k^^, > kj^^ > fc^ holds^under 
the experimentally reported load, r ~ 40pN-nm |T3|. The 
set of inequalities required between the chemical reaction 
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rates coincide with that for the pump process obtained 
here. Therefore the Fo part can work both as a motor 
and a pump efficiently by tuning the load r. In this sense 
these conditions are consistent. 

As far as we know, experiments which directly mea- 
sures the physical quantities, such as the rotation ve- 
locity and the pumping rate, have not been reported so 
far. However an experiment has been reported that in a 
whole FoFi ATP synthase the conformational change of 
one subunit of Fi part suppresses electronically and / or 
sterically the rotation hydrolysing ATPf^. This effect 
may be interpreted to vary the load effectively and con- 
trol which process occurs. 

In the model we considered here, only the diffusion 



process is taken into account. The highest efficiency ob- 
tained is e ~0.3, and it seems to be not high enough, since 
molecular motors are thought to transduce energy almost 
without loss There are several possibilities to solve 
this problem. The first is to reconsider the definition 
of the efficiency (E3j since our definition, following the 
argument in Ref.^^, is different from that in R.ef . [l^ . 
The second is to reconsider the relations between physi- 
cal quantities appropriate near the equilibrium, such as 
the Einstein's relation, since the system is far from equi- 
librium. The third is to introduce another mechanism, 
such as the electrostatic interaction between residues in 
Ref . [l5l| . These may lead to our better understanding of 
the energetics of non-equilibrium systems. 



[1] J.Howard, Mechanics of Motor Proteins and the Cy- 

toskelton (sinauer,2001) [10 
[2] R.D.Vale and R.A.Milligan, Science 288, 88 (2000). [11 
[3] R.D.Vale and F.Oosawa, Adv. Biophys. 26, 97 

(1990), A.Adjari and J.Prost, C. R. Acad. Sci. Paris [12 

II 315, 1635 (1992), N.J.Cordova, B.Ehmentrout, 

and G.F.Oster, Proc.Natl.Acad.Sci.U.S.A. 89 , 339 [13 

(1992), M.O.Magnasco, Phys.Rev.Lett. 71, 1477 (1993), 

R.D.Atsumian and M.Bier, Phys.Rev.Lett 72, 1766 [14 

(1994). 

[4] R.P.Feynman, R.Rayton, M. Sands, Feynman Lectures on [15 

Physics (Addison- Wesley, 1963). 
[5] F.Jiilicher, A.Ajdari, and J.Prost, Rev. Mod. Phys. 69, [16 

1269 (1997). 

[6] P.Reimann, Phys. Rep. 361, 57 (2002). [17 
[7] H.Linke, Appl. Phys. A. 75, 167 (2002) special issues on 

Brownian motors. [18 
[8] R.D.Atsumian and P.Hanggi, Physics Today 55(11), 33 

(2002). [19 
[9] S.B.Vik and B.J.Antonio, J.Biol.Chem. 269,30364 



(1994). 

Y.Sambongi et aL,Science 286, 1722 (1999). 
S.P.Tsunoda et al, Proc.Natl.Acad.Sci.U.S.A. 98, 898 
(2001). 

V.K.Rastogi and M.E.Garvin, Nature (London) 402, 263 
(1999). 

M.L.Hutchon et al, Proc.Natl.Acad.Sci.U.S.A 98, 8519 
(2001). 

W.Junge, H.Lill, and S.Engelbrecht, Trends. Biochem. 
Sci. 22, 420 (1997). 

T.Elston, H.Wang, and G.Oster, Nature(London) 391, 

510 (1998) and Supplementary Material. 

H.Miki, M.Sato and M.Kohmoto, Phys.Rev.E 68,061906 

(2003). 

S.P.Tsunoda et al, Proc.Natl.Acad.Sci.U.S.A. 98, 6560 
(2001) 

H.Noji, R.Yasuda, M.Yoshida, and K.Kinoshita, Na- 

ture(London) 386, 299 (1997) 

K.Sekimoto, J. Phys. Soc. Jpn. 66, 1234 (1997). 



